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ABSTRACT 

Recently, it was shown that the complex dynamical behaviour of spicules has to be interpreted as the 
result of simultaneous action of three kinds of motion: (1) field aligned flows, (2) swaying motions, and 
(3) torsional motions. We use high-quality observations from the CRisp Imaging SpectroPolarimeter 
at the Swedish 1-meter Solar Telescope to investigate signs of these different kinetic modes in spicules 
on the disk. Earlier, rapid blue-shifted excursions (RBEs), short-lived absorption features in the blue 
wing of chromospheric spectral lines, were identified as the disk counterpart of type II spicules. Here we 
report the existence of similar absorption features in the redwing of the Ca II 8542 and Ha lines: rapid 
red-shifted excursions (RREs). RREs are found over the whole solar disk and are located in the same 
regions as RBEs: in the vicinity of magnetic field concentrations. RREs have similar characteristics 
as RBEs: they have similar lengths, widths, lifetimes, and average Doppler velocity. The striking 
similarity of RREs to RBEs implies that RREs are a manifestation of the same physical phenomenon 
and that spicules harbour motions that can result in a net red-shift when observed on-disk. We find 
that RREs are less abundant than RBEs: the RRE/RBE detection count ratio is about 0.52 at disk 
center and 0.74 near the limb. We interpret the higher number of RBEs and the decreased imbalance 
towards the limb as an indication that field-aligned up-flows have a significant contribution to the 
net Dopplershift of the structure. Most RREs and RBEs are observed in isolation but we find many 
examples of parallel and touching RRE/RBE pairs which appear to be part of the same spicule. We 
interpret the existence of these RRE/RBE pairs and the observation that many RREs and RBEs have 
varying Dopplershift along their width as signs that torsional motion is an important characteristic 
of spicules. The fact that most RBEs and RREs are observed in isolation agrees with the idea that 
transversal swaying motion is another important kinetic mode. We find examples of transitions from 
RRE to RBE and vice versa. These transitions sometimes appear to propagate along the structure 
with speeds between 18 and 108 km s _1 and can be interpreted as the sign of a transverse (Alfvenic) 
wave. 
Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: corona 



1. INTRODUCTION 

The solar chromosphere is a highly dynamical atmo- 
spheric layer dominated by numerous small-scale struc- 
tures. At the solar limb, some of these small-scale dy- 
namical events can be seen as thin jet-like features pro- 
truding from above the photosphere to heights of 10 Mm 
or more. These jets are called spicul es and have b een 
known for more than a century, see Beckers (196 81) for 
a rev i ew of the older literature, and e.g., ITsiropoula 



et al. 



2012) for a more recent review. De Pontieu et al 



a| discovered the existence of a new class of spicules 
that differed from the known spicules by displaying only 
upwards motion, shorter lifetimes and higher ve l ocities . 



From an extensive set of data sets, |Pereira et al. ( 2012a[) 
determined typical velocities for these so-called type 11 
spicules to range between 30-110 km s _1 , and lifetimes 
between 50-150 s. The "other" type of spicules, or type 
I spicules, have typical lifetimes of 150-400 s, and show 
rising and falling phases in parabolic trajectories with 
maximum rise velocities ranging between 15-40 km s _1 . 
Type II spicules are the most abundant type which dom- 
inate both quiet Sun and coronal hole con ditions, only 
in active regions type I spicules dominate (De Pontieu 
et al.|2007a Pereira et al.|2012a). A recent investigation 



by Pereira et al. ( 2012b[ ) demonstrated that only modern 
observational techniques can provide the required spatio- 
temporal resolution that is necessary to detect these 
highly dynamical type II spicules and that bad seeing 
and technical limitations precluded earlier identification 
from the ground-based "classical" observations. 

Short lifetimes and high velocities are only two of sev- 
eral factors that contribute to make spicules notoriously 
difficult to observe: vigorous sideways and swaying mo- 
tion add to its dynamical complexity but the foremost 
facet that impedes unambiguous and continuous detec- 
tion throughout a spicule's lifetime is line-of-sight confu- 
sion. At and near the limb, the observer's line-of-sight 
crosses a large number of spicules that combine to a 
dense forest, in some spectral diagnostics almost a dif- 
fuse blanket, from which it is nearly impossible to trace 
all phases of a single spicule's evolution. For this rea- 
son, the identification of the disk counterpart of type 
II spicules was regarded as crucial for making progress 
on th e determ ination of its nature an d driving mecha- 
nism. i Langangen et al.| (12008 ) and Rouppe van der Voort 
et al. ( 2009J Paper 1) linked so-called "rapid blue-shifted 



excursions" (RBEs), short-lived and spatially elongated 
spectral asymmetries in the blue wings of the chromo- 
spheric Ca II 8542 and Ha lines, in on-disk observations 



to the type II spicules at t he lim b. Sekse et al. (2012 
Paper II) and Sekse et al. (2013, Paper 111) greatly ex- 
paneled on the statistics of these events and reinforced 
the interpretation of RBEs as the disk counterparts of 
type II spicules. 

Besides being one of the main constituents of the chro- 
mosphere, interest in type II spicules is motivated by 
their potential to play an important role in providing 
mass and energy to the upper solar atmosphere. Their 
characteristic sideways swaying motion has been inter- 
preted as a sign that the chromosphere is permeated 
with Alfvenic waves with sufficient energy to heat the 
quiet Sun corona and accelerat e the solar wind ( De Pon- 
tieu et al.||2007bD. Studies by |De Pontieu et al | ( |2U0^ T 
and De Pontieu et al. ( 2011[ ) provided evidence that type 
II spicules play an important role in mass loading and 
heating of the corona. The characteristic rapid fading 
observed in the Hinode Ca H observations has been in- 
terpreted as a sign of heating of the structure to tempera- 
tures that leaves no significant opacity in the observation 
passband. 

As mentioned above, type II spicules display a complex 
dynamical behaviour. Besides apparent rising velocities 
that vary over a range of more than 100 km s _1 , many 
are displaced transversally over distances up to 1 Mm 
and may reach trans verse velocities up to 60 km s _1 
(Pereira et al. 2012a). Oscillatory swaying motion has 



been observed with amplitud es of order 10-20 km s 1 
and periodicities of 100-500 s ([ De Pontieu et al |2007b[ ). 
Okamoto & De Pontieu (2011) found a mix ol upward 
and downward propagating as well as standing waves in 
a large number of spicules in high-cadence Hinode Ca H 
limb observations. Similar periods and amplit udes were 
found for Ho ; RBEs at disk center in Paper III. Suematsu 
( |2008| ) reported that some spicules appear as do 
d sometimes even multiple threads with evidence 
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spinning motion. Judge et al. (2012) reported the strik- 
ing observation of spicule- like features that suddenly ap- 
pear along their full length in extremely high-cadence Ha 
wing observations close to the limb. Paper III reported 
similar behaviour for RBEs at disk center. If upward 
moving plasma were the only motion in spicules, sudden 
appearance would imply extreme propagati on speeds in 
excess of several thousand km s _1 . However, De Pontieu 



et al. (2012) established that besides field-aligned flows 
and swaying motions, a third type of motion, torsional 
motions of order 25-30 km s _1 are governing the dynam- 
ics of spicules. When torsional motions are considered, 
the sudden appearance of spicule-like features in fixed- 
wavelength imaging can be explained as the result of ro- 
tational Dopplershift into the observed passband (Pa- 
per III). The observed omni-presence of torsional motion 
in spicules implies that the combined action of all three 
modes of motion needs to be considered when interpret- 
ing Doppler signals to estimate the transport of plasma 
and Alfvenic wave energy by spicules through the solar 
atmosphere. The ubiquity of torsional motion may also 
have important implications for determining the produc- 
tion mechanism of spicules and suggest that they may be 
an important agent for transporting helicity to the upper 

solar atmosphere. 

In th is study, we continue the work by |De Pontieu et aT 
(2012J) and Papers I— III to investigate signs of torsional 



motion s in the disk coun t erpar t of type II spicules. The 
work of De Pontieu et al. (2012 ) was based on various slit 
and imaging spectroscopic data sets at the limb. Juxta- 
position of the three kinetic modes was observationally 
established by curved spectral profiles in off-limb Ca II H 
X-t diagrams, and tilted spectral profiles in X-x diagrams 
of the Ca II H, Ha, and Ca II 8542 spectral lines. In 
addition, the striking similarity of rows of on-disk and 
near-limb absorption features in red and blue wing sam- 
plings was noted as an important indication of the pres- 
ence of large-amplitude transversal velocities. In the blue 
wing, these absorption features are comparable to RBEs 
observed closer to disk center. Observed in the red wing, 
these absorption features are more appropriately referred 
to as "rapid red-shifted excursions", or RREs. Here we 
report the presence of RREs over the whole solar disk 
and characterise their properties. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were obtained at the Swedish 1-m So- 
lar Telescope (SST, |Scharmer et al.||2003a[ ) on La Palma 
using the CRisp Ima ging SpectroPolarimeter (CRISP, 



Scharmer et al. 2008) instrument. The CRISP instru- 
ment is based on a dual Fabry-Perot interferometer (FPI) 
and contains three high-speed low-noise CCD cameras 
which operate at a frame rate of 35 frames per second 
with an exposure time of 17 ms per frame. All three 
cameras are located behind the CRISP pre-filter and are 
synchronised by an optical chopper. Two of the cam- 
eras are also positioned behind the FPI and a polarising 
beam splitter, while the third camera, which acts as an 
anchor channel for image processing, is located before 
the FPI. CRISP has a field of view (FOV) of approxi- 
mately 61 // x61 // and an image scale of / /0592/pixel. Fast 
wavelength tuning (<50 ms) is possible between any two 
positions within the spectral region given by the spectral 
width of the pre-filter. With its high speed capabilities 
CRISP is an ideal instrument for studying the dynam- 
ics of the chromosphere through imaging spectroscopy. 
Here we are interested in the Ca II 8542 and Ha spec- 
tral lines. For Ca II 8542, CRISP has a transmission 
FWHM of 111 mA and a pre-filter FWHM of 9.3 A. 
The transmission FWHM for Ha is 66 mA while the pre- 
filter has a FWHM of 4.9 A. The combination of the SST 



adaptive optics system (Scharmer et al. 2003b) and the 
Multi-Object, Multi-F rame Blind Deconvolut ion image 
restoration technique ( |van Noort et al."|2005[ ) results in 
high spatial resolution, down to the telescope diffraction 
limit (X/D=0'!l4 for Ha). Remaining small-scale see- 
ing deformations introduced by the non-simultaneity of 
the narrowband CRISP images is mini mized by using the 
cross-correlation method developed by He nriques] (|2012 ) . 
We analyze several high-quality datasets, each with 
different properties so that various properties of RREs 
can be addressed. Table [I] provides details on the dif- 
ferent data sets. Data sets A and B are time series 
with extended wavelength coverage of both wings of the 
Ca II 8542 and Ha spectral lines. The other time series of 
2011 May 5 (data set C) has an extremely high cadence 
(0.88 s) so that the dynamical evolution of RREs and 
RBEs are temporally resolved. Only 4 wavelength po- 
sitions are sampled, ±1032 mA, line center, and the far 
blue wing at —2064 mA. For more details on the 2011 



Table 1 

Observational data sets. 



Data set 


Date 


Spectral line 


# Line positions 


Cadence 

(s) 


Duration 
(h:m:s) 


Target 


/x 


A 


2008 Jun 15 


Ca II 8542 


28 


11.33 


00:40:25 


CH 


0.987 


B 


2011 May 05 


Ra 


35 


7.96 


00:52:31 


QS 


1 


C 


2011 May 05 


Ra 


4 


0.88 


00:30:53 


QS 


1 


D 


2011 Sep 18 


Ca II 8542 


47 


10.37 


02:18:53 


QS 


1 


E 


2012 May 26 


Ra 


33 


(...) 


(...) 


QS, CH 


[1-0] 


F 


2012 Jun 02 


Ra 


33 


(...) 


(...) 


QS, CH 


[1-0] 



Note. — CH: coronal hole, QS: quiet Sun. 



May 5 data sets, see Paper III. Data set A was anal- 
ysed in Paper I but the data was re-processed for this 
study: now an efficient correction was applied to deal 
with artefacts caused by the semi-transparency of the 
back-illum inated CCD chips at lo nger wavelengths (for 
details see de la Cruz Rodriguez 2Q1Q[ ). Artefacts were 
effectively removed and led to a significant increase in 
the number of RBE detections. 

For the other 3 data sets, the observations were 
recorded under conditions of "snapshot" seeing: the 
seeing varied between extremes of total blurriness and 
crystal-clear sharpness. Most of the time during snap- 
shot seeing, the conditions are useless for science anal- 
ysis but the intervals of good quality are of sufficiently 
long duration to cover the acquisition of a complete spec- 
tral line scan. On 2011 September 18 (data set D), this 
yielded a large number of disk center Ca II 8542 line scans 
with a dense sampling down to 55 mA (half of the CRISP 
transmission FWHM). On 2012 May 26, Ha line scans at 
7 disk positions were acquired (data set E), starting at 
disk center and following the central meridian down to 
the South Pole in intervals of \i = 0.1 (with \i = cos(9) 
and the observing angle). On 2012 June 2 (data set F), 
the central meridian from disk center to the North Pole 
was covered by 22 Ha line scans (down to 86 mA sam- 
pling), combining into a near-continuous center-to-limb 
scan (continuous with the exception of one narrow inter- 
ruption - see Fig [7]). Precise pointing information was 
obtained from co-alignment with full disk images from 
the 1 700 A channel of the AIA instrument (Lemen et al. 
2012[ ) on NASA's Solar Dynamics Observatory. Using 
photospheric bright points in the far Ha wing as refer- 
ence, sub-arcsecond accuracy in the co-alignment can be 
achieved. 

3. METHOD 

For each dataset used in this study, an automated de- 
tection routine designed for locating and analyzing RBEs 
was used to find both RBEs and RREs. The same detec- 
tion metho d (with slight modifications) was used in ear- 
lier studies ( Rouppe van der Voort et al. 2009 Sekse et al 



2012[ |2013| ) When searching lor KKEs, the automated 
detection routine worked on inverse Doppler images and 
hence, found RREs in a similar way as RBEs. For Ha, 
the Dopplergrams were at ~45 km s _1 , for Ca II 8542 
at 20 km s _1 . In the process of detecting RREs we set 
the threshold of minimum length of the detected events 
to 20 pixels or ~0.86 Mm (17 pixels for data set A from 
2008, when the images scale was slightly larger then for 
the later data sets). This is longer than for our previ- 



ous studies (17 pixels in Paper II and 14 pixels in Pa- 
per III) which was found to be necessary in order to 
avoid detecting false positives. For example, false pos- 
itives are caused by the presence of small-scale down- 
flowi ng features, or flocculent flows (Vissers & Rouppe 
|van der Voort||2012| ). Flocculent flows are ubiquitously 
present in the vicinity of magnetic field concentrations 
and have distinctively different morphology and dynami- 
cal behavior than the RREs that we study here. For each 
detection, the automated routine records the spatial and 
temporal position as well as determining the Doppler ve- 
locity and Doppler width along the event. Doppler ve- 
locity and width were determined on a residual profile, 
obtained from subtracting the RRE/RBE profile from 
a reference profile that was constructed from averaging 
line profiles over the whole FOV. Furthermore, the au- 
tomated detection routine attempts to link up detected 
events in successive frames into multi-frame RREs and 
RBEs which can be used for measurements of lifetimes, 
transverse motion, and longitudinal motions. 

Every dataset was analysed for both RREs and RBEs 
in order to make comparisons and previously searched 
datasets were re-examined with the updated length con- 
straint in order to make a proper comparison. 

The verification of RREs and the creation of space- 
time diagrams was done with the widget based anal- 



ysis t ool CRISPEX ([Vissers & Rouppe van der Voort 



2012[ ) which allows for efficient exploration ol multi- 
dimensional datasets. 

4. RESULTS 

When searching on-disk Ha and Ca II 8542 datasets 
for excursions in the red wing, many events are detected. 
These events are very similar in appearance and mor- 
phology as RBEs. Figures [T| and [2] show details of RREs 
in comparison with RBEs. Figure]!] shows an RRE event 
observed in Ha. Just after the event disappeared in the 
Ha red wing, an RBE appeared in the blue wing just 
adjacent to where the RRE was before. This RBE is 
shown in the right column of Fig. [I] Figure [2] shows 
RRE and RBE events in Ca II 8542. These events were 
visible as an RRE/RBE pair of adjacent events of sim- 
ilar length and morphology giving a strong impression 
of being part of the same spicule structure. The mid- 
dle panels (c and d) of these figures show the detailed 
line profiles at a central location in the event, together 
with reference average line profiles. Just like for RBEs, 
the RREs are characterised by an asymmetry of the line 
profile with enhanced depression of the line wing - in 
the red wing rather than the blue wing. Likewise, the 
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Figure 1. Details of an RRE (left column) and an RBE (right 
column) that occur in rapid succession in neighbouring locations 
in Ha (data set B). The top row shows the temporal evolution of 
the spectrum at a central location in the RRE and RBE, marked 
with asterisks in panels (e) and (f). The middle row shows the 
detailed spectral profiles (solid line) of the positions marked with 
asterisks in panel (e) and (f ) at the time of the background image. 
The line profile marked with diamonds is an average over the whole 
FOV shown as reference. The thick dashed line is the difference 
between the average profile and the RBE/RRE profile which is used 
to calculate the Doppler velocity and width of the event from. The 
bottom row shows 35 km s — 1 Doppler images of the RRE (dark) 
and RBE (bright) separated by a couple of timesteps. 



temporal evolution for RREs and RBEs is similar: the 
asymmetry of the line profile is rather short-lived with 
comparable lifetimes for the two type of events. This is 
illustrated in the top panels (a and b). Furthermore, as 
can be seen from the bottom panels, the morphology for 
RREs and RBEs is similar: narrow, elongated structures 
often protruding from regions with enhanced magnetic 
fields (bright points, network regions, etc). 

The measurements of Doppler velocities, Doppler 
widths, and lengths for detected events in the data sets A 
(Ca II 8542) and B (Ha) are shown in Fig. [3] As a way to 
suppress the contribution from possible false detections, 
we only include measurements on events that are part 
of a multi- frame event, i.e., events that could be iden- 
tified in successive detection images and were therefore 
connected over a sequence of time steps. These mea- 
surements were done per time step. Doppler velocities 
and width measurements are averages over the detection 
length of the event. There are significantly more events 
in the blue wing than in the red wing: for the 52 min Ha 
quiet Sun time series (B) there are 3128 RBE and 1354 
RRE detections (43%), whereas for the 40 min Ca II 8542 
coronal hole time series (A) there are 2544 RBE and 1237 
RRE detections. 

From the histograms in Fig. [3] it is clear that for these 
properties, RREs (solid lines) are very similar to RBEs. 
Each property shows largely overlapping distributions. 
For Ca II 8542 RREs we see Doppler velocities lying 




Figure 2. Details of an RRE and RBE observed in Ca II 8542 
(data set A). The panels follow the same format as Fig. fl] Both 
the RRE and RBE are present in the same timestep and appear to 
be part of the same spicule structure. The background 20 km s _1 
Dopplermaps in panels (e) and (f) are the same. The same event 
is shown in more detail in the right column of Fig. [8] 



mainly between 10 and 20 km s _1 , Doppler widths rang- 
ing from 1 to 20 km s _1 , and lengths ranging mainly 
from 1 Mm up to approximately 4.5 Mm. These RRE 
properties are comparable to the RBE properties (dashed 
lines) found in the same data set with only minor differ- 
ences in that Doppler widths and lengths for RBEs are 
shifted to slightly higher values: average Doppler width 
for RREs at 7.0 km s" 1 and for RBEs at 8.2 km s _1 , 
and average length for RREs at 2.5 Mm and for RBEs 
at 2.9 Mm. For Ha, the distributions for Doppler veloci- 
ties (ranging between ~20 and ~35 km s _1 ) and lengths 
(between ~1.5 and ~4.5 Mm) are almost identical for 
RREs and RBEs. For Ha Doppler widths, the distribu- 
tion for RREs is slightly lower: the average for RREs is 
7.5 km s _1 and for RBEs 8.7 km s" 1 . 

The variation of Doppler velocity and Doppler width 
along RREs and RBEs is shown in Figs. [I] (Ha) and [5] 
(Ca II 8542). In these figures, the variation of Doppler 
signal is shown in colour coding drawn at the physical 
locations of RREs and RBEs on a background image of 
the FOV. For each pixel, colour indicates the average 
value for all events that occurred at that location during 
the time series. Comparing the top and bottom rows 
in these figures, it is evident that there are more events 
in the blue wing than in the red wing: RBEs are more 
frequent than RREs (see above for detailed numbers). 
Furthermore, the spatial distribution of RREs is very 
similar to RBEs: mostly concentrated around regions of 
enhanced magnetic activity. RREs and RBEs occur at 
the same spatial locations. 

Paper III shows comparable maps of Ha RBE Doppler 
velocity and width for the same data set (B) in their 
figure 1 as we do here in Fig. ^ and d. The differ- 
ences are that in Paper III a shorter detection length 
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Figure 3. Histograms of RRE (solid line) and RBE (dashed line) properties in Ca II 8542 (top row) and Ha (bottom row). Doppler 
velocities (a and d) and Doppler widths (b and e) are averaged over the length of the event. Lengths are presented in panels (c) and (f). 
Number of measurements, Ca II 8542 data set A: 1237 RREs and 2544 RBEs, Ha data set B: 1354 RREs and 3128 RBEs. 



constraint was used (resulting in 1.7 times more RBE de- 
tections) and a slightly different method of drawing these 
maps: events were drawn in sequential order so that later 
events mask earlier events. The maps in Fig. ^p and d 
have smoother appearance than their counterparts in Pa- 
per III but the same picture emerges: there appears to 
be no clear general trend for the variation of Doppler 
properties along the length of RBEs in this disk center 
quiet Sun time series. Similarly, the RRE Doppler prop- 
erties in the top row panels seem to lack a clear general 
trend too. In the Doppler velocity map of Fig. [4^i, there 
might be a weak preference for higher velocities towards 
the footpoints, closer to the magnetic network concen- 
trations (a trend of more red in the lower parts of the 
RREs). For the Doppler widths in Fig. Eb, there might 
be a preference for lower widths towards the footpoints 
(more blue towards the bottom). Note that the trend of 
more dark and blue colours in the RRE Doppler width 
map as compared to the RBE Doppler width map re- 
flects the lower Doppler width distribution for RREs in 
Fig#. 

The Doppler property maps for Ca II 8542 RBEs 
in Fig. ^jp and d can be compared with the left pan - 
els of figure 12 of Rouppe van der Voort et al. (12009). 
They showed Doppler properties of 413 Ca 11 8542 RBEs 
from the same data set but detected in Dopplermaps at 
30 km s -1 . Here we show results from 3655 Ca II 8542 



RBE events detected in 20 km s _1 Dopplermaps. We 
confirm the results from Paper I: in this coronal hole 
data set, there is a trend for increasing Doppler velocity 
and width along RBEs (in the colour scheme of this fig- 
ure: a trend of more red towards the top end of RBEs). 
For Ca II 8542 RREs in the top panels of Fig. [5] the vari- 
ation is much more diverse and it is hard to see any clear 
trend. 

4.1. Center to Limb Variation 

To further explore the occurrence of RREs on the solar 
disk, we ran the detection method on data sets at various 
disk positions. The detection was done on Dopplergrams 
at about the same wavelength offset for all different data 
sets: approximately ±45 km s _1 . Figure [6] shows the 
number of RRE and RBE detections as function of disk 
position p = r /R@ with r the distance from disk center 
and i?0 the solar radius. The numerous dots around 
p = are the detection numbers from every frame of 
the 52 min data set A taken at disk center (396 time 
steps). Generally, more RBEs are detected than RREs. 
Although there is a significant scatter to the number of 
detections, we see a clear trend of increasing numbers of 
both RREs and RBEs being detected towards the limb. 
The solid and dashed lines are linear fits to the RRE 
and RBE detection numbers, respectively, and show an 
increase in the average number of detections from 5.3, at 
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Figure 4. Ha Doppler properties as function of position along RREs and RBEs drawn at their location in the FOV. RREs are plotted 
in the top row and RBEs in the bottom row, with Doppler velocities being displayed in the left column and Doppler widths in the right 
column. Background images are ±1032 ml Doppler image (a), line core (b), +1032 mA red wing (c), and far blue wing image (d). 



disk center, to 14.4, near the limb, for RREs, and from 
10.1 to 19.4 for RBEs. 

Figure [7] shows a number of close-up images centered 
on network regions at various disk positions. The center- 
to-limb scan of data set F mainly covered quiet Sun, it 
crossed a large equatorial coronal hole between heliocen- 
tric y=0 // -210 // . The upper-left close-up image is from 
the coronal hole (y=152 // -174 // ). Over the whole disk, 
RREs and RBEs are observed in close vicinity. Most 
often, RREs and RBEs are seen in isolation, but there 
are a number of examples of RRE/RBE pairs where the 



RRE and RBE signal is observed in parallel and touching 
structures, suggestive of being part of the same physical 
entity. A prominent example is visible in the lower left 
panel (?/=694 // -704 // ) but close inspection reveals exam- 
ples in every close-up image. 

4.2. Torsional motion 

Parallel and touching RRE /RB E pairs are found over 
the whole solar disk (see Sect. 4.1 ). Figure M shows exam- 
ples from high spectral resolution Ca II 8M2 scans close 
to disk center. The middle panels show spectrograms (A- 




Figure 5. Ca II 8542 Doppler properties as function of position along RREs and RBEs drawn at their location in the FOV. RREs are 
plotted in the top row and RBEs in the bottom row, with Doppler velocities being displayed in the left column and Doppler widths in the 
right column. Background images are Ca II 8542 ±582 mA Doppler image (a), line core (b), +582 mA red wing (c), and far blue wing (d). 



x diagrams) from an artificial "spectrograph slit" drawn 
in the direction transversal to the main axes of selected 
events close to a network region. The asymmetry of the 
spectral profiles changes rapidly from one wing to the 
other over short distances, resulting in tilted profiles in 
the spectrograms. For the main event in the left column, 
centered on 0" along the slit in the middle panel, the 
asymmetry changes from approximately +20 km s _1 to 
—20 km s _1 over only 250 km. Another less extreme 
example is visible at — 3.5" along the slit. The right col- 
umn shows an example of an RRE/RBE pair with the 



asymmetry varying between approximately ±20 km s 1 
over the width of the structure. 

4.3. Temporal change between RREs and RBEs 

Besides parallel and touching RRE/RBE pairs, we fre- 
quently observe events for which the asymmetry changes 
from red-ward to blue-ward and vice versa. These 
changes are very rapid and the most convincing exam- 
ples that these modulations really are occurring in the 
same physical structure can be found in the high tempo- 
ral resolution dataset C (cadence 0.88 s). Figure [9] shows 
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Figure 6. Number of detections of RREs (red) and RBEs (blue) 
per time step over the FOV of several datasets going from disc 
center (p = t/Rq = 0) to the solar limb (p = 1). Corresponding 
p = cos(observing angle) scaling is shown on top. Dots near p = 
signify data points from the disk center data set B, the spread in p 
is artificial to show the density of data points. Stars are from data 
set E, crosses are from F. The histogram bars are averages over the 
data points in Ap = 0.1 intervals. Linear fits to the data points 
are shown with a solid (RREs) and dashed (RBEs) line. 



a few examples. The top panels show an event that starts 
as a long RRE of about 3" length which changes into a 
short RBE that appears to be growing from the bottom. 
While the RBE appears, the RRE vanishes. This change 
occurs at a propagation speed that is estimated to be 
about 108 km s -1 . The bottom panels show two events 
that appear close to each other. The first event starts as 
an RRE (dark feature) that is centered on ~ {11'! 5, 34'/ 5) 
at t — —16.1 s and moves slowly (diagonally) upwards. 
At the same time, an RBE slowly grows from the bottom 
with a slight acceleration after about t = 5 s. The av- 
erage propagation speed of the change is approximately 
18 km s -1 . Almost parallel to this event, there is a simi- 
lar event that starts as an RRE that is not visible before 
t = s, the space-time diagram is shown as the right- 
most panel. While the RRE slowly moves upward, an 
RBE grows at an approximate speed of 48 km s -1 . Sim- 
ilar space-time diagrams as the right panels of Fig. [9] were 
made for 7 other examples, finding propagation speeds 
between 18 and 102 km s _1 . The average over all events 
is 58 km s _1 . 

5. DISCUSSION 

The red wing counterpart of RBEs: narrow, elongated 
and short-lived absorption features in the Ca II 8542 and 
Ha lines in the form of asymmetries in the red wing, 
or RREs, can be observed in large abundance over the 
whole solar disk. They are very similar to RBEs in terms 
of length and average Doppler velocity, while the average 
Doppler width is slightly lower for RREs than for RBEs. 
Their spatial distribution is similar to RBEs: they are 
found around concentrations of enhanced magnetic field. 
We did not do a detailed study of the temporal evolu- 
tion of RREs like in Paper III, but cursory inspection of 
our data suggests that the lifetimes of RREs are similar 
to RBEs. RREs are less abundant than RBEs, around 
disk center about half the number density, but significant 
numbers can be found anywhere on the disk. 

When we reflect on our own work on the disk counter- 



parts of type II spicules, we can identify several reasons 
why RREs were not identified earlier. One reason is that 
the observations from the early CRISP campaigns were 
not optimally suited to identify RREs. Limiting factors 
for the two data sets from 2008 that were used in Paper I: 
the Ha time series had an asymmetric spectral sampling 
that did not cover the red wing further than 800 mA. The 
Ca II 8542 time series (data set A in this study) suffered 
from artefacts related to the sem i-transparency of the 
CCD chip at long wavelengths (see |de la Cruz Rodriguez] 
2010). For this data set, this problem was not resolved 
until 2011. The 2010 observations acquired for follow-up 
studies (used in Paper II) were specifically designed for 
RBE detection: the Ha line was only sampled in the blue 
wing, and the Ca II 8542 sampling was relatively coarse. 
It was not until high resolution off-limb Ca II H slit 
spectroscopy could be combined with high quality CRISP 
limb observations that the ubiquitous presence of tor- 
sional motions in spicules could be unequivocally estab- 



lished ( |De Pontieu et al. |2012| ). It was realised that 
spicules undergo three different types of motions at the 
same time: (1) field- aligned flows of order 50-100 km s _1 , 
(2) swaying motions of order 15-20 km s _1 , and (3) tor- 
sional motions of order 25-30 km s _1 . The presence of 
these three kinetic modes explained the appearance of 
spicular features in near-limb blue and red wing CRISP 
Ha and Ca II 8542 images. Bushes of red-shifted and 
blue-shifted features appear with similar morphology and 
orientation. Field-aligned flows towards and away from 
the observer alone cannot explain the striking similarity 
as the asymmetry of near-limb line-of-sight projection 
should cause large morphological inequality. Doppler 
modulation due to transverse motions from modes (2) 
and (3) can explain the presence of near-cospatial and 
morphologically similar absorption features in the blue 
and red wings of the near-limb CRISP observations. This 
understanding led us to this work investigating the pres- 
ence of similar red-wing absorption features at disk po- 
sitions further away from the limb. 

We interpret the presence of RREs close to disk cen- 
ter as a strong indication that field-aligned flow cannot 
be the only type of plasma motion in spicules. High- 
amplitude transversal motions must be present too in 
order to produce an effective Dopplershift into the red 
wing. We find evidence for these transversal motions to 
be in the form of torsional motions and/or in the form 
of swaying motions. 

The observation of parallel and touching RRE/RBE 
pairs is a direct indication of torsional or spinning motion 
in spicules (see Figs.[2j [6l and|8|. We find a significant 
number of RRE/RBEpairs in our data but most com- 
mon are isolated RREs and RBEs. The spectral profiles 
in X-x diagrams extracted perpendicular to the main axes 
of isolated events in high-spectral resolution observations 
often show signs of tilts similar (but less pronounced) to 
the extreme cases of the RRE/RBE pairs in the mid- 
dle row of panels in Fig. [8j This is most clear in the 
Ca II 8542 line which has a narrower core than Ha and 
therefore more pronounced Doppler modulation. This 
suggests that torsional motion is likely to be present in 
a large fraction of spicules and therefore an important 
kinetic mode. 

The importance of swaying as a kinetic mode in RBEs 
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Figure 7. Examples of RRE ± RBE concentrations around network regions at different disk positions from the center-to-limb scan of data 
set F. The left panel shows the various telescope pointings as white squares on the background of a co-temporal SDO/AIA 304 A image. 
The 9 images are cutouts from ±1032 Ha Dopplergrams, each centered on network concentrations with enhanced RRE/RBE activity. 
RREs are white, RBEs are black. Axis labels in heliocentric coordinates (arcseconds). 
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Figure 8. Examples of direct evidence of torsional motion in the 
disk counterpart of spicules. Spectral profiles are extracted from 
Ca II 8542 datasets (data set D in left column, data set A in the 
right column) in a direction perpendicular to the main axes of 
a number of RREs and RBEs (the "slit" drawn as white dashed 
lines in the Ca II 8542 images in the top panels). Top panels show 
Ca II 8542 red wing images (left: +440 mA, right: +582 mA). 
Corresponding spectrograms (or X-x diagrams) are shown in the 
middle panels, with the origin of the y-axis centered on the main 
RRE/RBE pair. The bottom panels show ±16 km s _1 (left) and 
±20 km s _1 (right) Dopplergrams with RREs black and RBEs 
white. Figure [2] shows more details of the right hand-side event. 



was already evident from transverse motion measure- 
ments (see Papers I— III) . The observation that most 
events are either RRE or RBE and not pairs can be in- 
terpreted as swaying producing significant offsets in the 
Dopplershift of these structures. Field- aligned up-flows 



can also result in a Doppler offset so that the dominance 
of isolated RRE and RBE events is probably due to the 
combined effect of both kinetic modes. Figure 10 illus- 
trates the effect of swaying and/or up-flows on the ob- 
served Dopplershift from a tube with torsional motion. 
Only in the situation that these modes do not produce a 
significant net velocity component along the line-of-sight, 
a clear RRE/RBE pair can be detected. 

The observation of finding generally more RBEs than 
RREs suggests that plasma up-flows are present in 
spicules and result in a bias towards blue-shifts. For 
the center-to-limb variation of the detection numbers, 
we see that the total number of detected events increases 
towards the limb (i.e., both more detected RBEs and 
RREs). In addition, we see that the number of RREs in- 
creases relatively more than the number of RBEs. Partly, 
the increase in total number of detected events towards 
the limb can be explained by the curvature of the Sun 
and the related increase of the solar surface covered by 
the fixed size of the FOV when pointing closer to the 
limb. In addition, this can be interpreted as resulting 
from transversal motions gaining weight as the line-of- 
sight becomes more inclined to the average spicule axis 
as we go towards the limb. The net contribution from 
field-aligned flows would decrease closer to the limb and 
we see a more balanced ratio of RRE to RBE detections: 
from RRE/RBE « 0.52 near disk center to « 0.74 near 
the limb. 

The transition from RREs into RBEs and vice versa 
(see Fig. [9| can be interpreted as the result of an Alfvenic 
wave associated with the swaying motion. The exact na- 
ture of the wave is difficult to determine since the asso- 
ciated time scale is so short (high propagation speed). 
Only in the high temporal resolution dataset (cadence 
0.88 s) we found a number of convincing examples where 
these changes were occurring in the same physical struc- 
ture. Such high cadence observations could only be ac- 
quired at the expense of high spectral resolution so that 
valuable diagnostic information that could provide more 
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Figure 9. Examples of the rapid change from RRE to RBE in the fast Ha dataset C. Left panels: sequential images showing the time 
evolution of the transition from one form into the other. The right panels show the space-time (x-t) diagrams from a slit drawn along the 
RRE/RBE axis. The dashed lines serve as reference for a propagation speed for the transition. In the bottom panels, two examples of 
transitions from RRE to RBE can be seen. Both events are diagonally oriented in the images on the left. The left x-t diagram is associated 
with the RRE (dark feature) that is visible from the start of the sequence of images. The right x-t diagram is associated with the event 
slightly below that does not appear as an RRE before t = 0. 

insight into the physical nature of the changes is lacking. the solar atmosphere. 



6. CONCLUSION 

In a number of high-quality data sets with different ob- 
servational characteristics (variation in e.g., spectral line 
sampled, temporal and spectral resolution, observing an- 
gle and target) we find clear evidence of the presence of 
three different types of motion in the disk counterpart 
of type II spicules: (1) field aligned flows, (2) swaying 
motions, and ( 3) torsional motions. Ou r results support 
the findings of De Pontieu et al. (2012) who established 
the presence and relative importance of these three ki- 
netic modes in spicules at the limb. Rotational motion in 
spicules had been reported in the earlier lit erature (e.g., 
Beckers|1968l|Pasachoff et al.|1968||Weart|1970| ) but only 
with the aid of modern observational techniques the ubiq- 
uity of this type of motion could be determined. We also 
see evidence for propagation of both the swaying and 
torsional motions, which is compatible with the presence 
and propagation of Alfvenic waves. Our analysis on the 
disk cou nterpart of type II s picule s reinforces the conclu- 
sions of De Pontieu et al. (2012): the combined action 
of all three kinds ot motion needs to be considered when 
interpreting spicule Doppler signals to estimate their con- 
tribution to the transport of plasma and energy through 
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